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Abstract: It is shown that for any elastic string model with energy density U and tension 
T, the divergent contribution from gravitational self interaction can be allowed for by an 
action renormalisation proportional to {U — T)'^. This formula is applied to the important 
special case of a bare model of the transonic type ( characterised by a constant value of the 
product UT ) that represents the macroscopically averaged effect of shortwavelength wiggles 
on an underlying microscopic model of the Nambu-Goto type (characterised by U = T). 



1 Introduction 

Although not so important for hghtweight cosmic strings, such as may may have been 
formed later on (for example at the time of electroweak symmetry breaking) gravitational 
self interaction is generally supposed to have played an essential role in the evolution of 
the cosmic strings that may have been formed earlier, at the epoch of G.U.T. symmetry 
breaking. In the kind of scenario Q originally proposed by Kibble an initial period during 
which the main damping mechanism was the friction exerted by the ambient thermal gas 
would be followed by a period during which the main damping mechanism would have 
been gravitational, at least for the local kind of strings to be considered here (in contrast 
with global strings for which axion radiation damping would have been more important). 

Nearly all the work that has been done Q on gravitational self interaction in cosmic 
strings has been based on the use of a string model of the simplest kind, namely that of 
Nambu-Goto, in which (with the speed of light set to one) the energy density U and the 
tension T are both equal to a constant m^, where m is a mass scale that will typically be 
of the order of the mass associated with the Higgs field responsible for the relevant vacuum 
symmetry breaking. The use of such a simple model will be justifiable while gravitational 
self interaction is important, even in typical cases where the strings are of the current 
carrying kind whose likely relevance was first pointed out by Witten since, at least at 
the outset, the currents would be expected to be very weak so that their effects would be 
relatively negligible. 
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At a later stage the currents in small loops would intensify as the loops contracted due 
to radiative energy loss, so that it might become necessary to use a description based on an 
elastic string model of an appropriate kind j^, in terms of which it is possible to describe 
oscillations [Q, ||] about stationary "vorton" states (whose existence would not even 
be possible if the Nambu-Goto description remained valid). However by this later stage 
(except in the case of an electromagnetically neutral current) the main self interaction 
mechanism would no longer be gravitational but electromagnetic. Thus so long as one has 
to do with a regime in which gravitation is still the main self interaction mechanism, a 
Nambu-Goto description will usually be adequate. 

Although adequate in principle, an absolutely precise Nambu-Goto description will 
seldom be feasible in numerical simulations, due to the enormous amount of information 
that would be required for following the detailed evolution of small scale wiggles. Some 
kind of approximation will therefore be needed in pratice. Following the remark by Shellard 
and Allen 0] that in the presence of small scale wiggles the macroscopically averaged value 
of the energy density U would exceed its Nambu-Goto value m^, I observed § that the 
tension T would be correspondingly diminished, and proposed as an approximation the 
use of a model in which the product of energy density and tension remained constant, 

UT = m^. (1) 

This relation was confirmed for a special subclass of wiggle modes by Vilenkin |^ and more 
recently for unrestricted wiggle modes by Martin ||l^ . A particularly attractive feature of 
this model is the property of being transonic, in the sense that the propagation speed 

Ce = \/t/U (2) 
of its extrinsic (tranverse) perturbations is the same as the speed 

c, = ^-dT/dU (3) 

of its longitudinal (sound type) modes. This makes it possible to prove that such a 
model can match the long term evolution of the string with any desired acuracy (depending 
on the resolution) and without any accumulative error buildup, at least in flat spacetime 
where the dynamical equations are exactly integrable. 

For the practical description of cosmic strings in circumstances where the main kind 
of self interaction is gravitational it is therefore this transonic model that will commonly 
be most appropriate. The purpose of this article is to consider the way this model will 
need to be modified to allow for the dominant effect of the self interaction, which will be 
divergent. 

As in the simpler case of electromagnetic self interaction, gravitational self interaction 
in a four dimensional spacetime background will give rise in point particle models to 
pole type singularities, and in string models to logarithmic singularities, which need to be 
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regularised by the use of an "ultraviolet" cut off. In order to do this in a systematic manner 
I have developed a simple technical formalism |12| whose application to the electromagnetic 
self interaction gives results that agree with the conclusions of earlier approaches, while 
providing a more convenient means of dealing with applications to particular problems 
such as oscillations of a conducting string loop about a stationary vorton state p3[|. The 
more recent application of this formalism to the linearised gravitational case |jlj] has 
however given results [15| that deviate from what had previously been obtained in 
the only case that had been considered previously, namely that of a Nambu-Goto string, 
for which the divergent part simply vanishes. A detailed examination of this particular 
case 1 16] has shown that previous assertions to the contrary Q were effectively based on 
the neglect of terms that would indeed be relatively small when due to high frequency 
gravitational radiation from a distant source, but that in the case of local self gravitation 
are of the same order as the other terms, which they finally cancel. 

What will be done here is to apply the general formula] 15 1 for the divergent part of the 
self gravitational interaction of a string to the non trivial case of a generic elastic string 
model, and in particular to the case of the transonic string model, for which it does not 
vanish. The non vanishing result for the transonic string model is entirely consistent with 
its interpretation as a smoothed approximation to a wiggly Nambu-Goto model for which 
the corresponding short range self gravitational contribution vanishes. The non-vanishing 
short range contribution in the approximate description represents the finite intermediate 
range contribution in the underlying Nambu-Goto model. Thus the treatment provided 
here takes care of everything except the very long range part of the self interaction, which 
will be dynamically negligible in the short run, though it is of course very important in 
the long run since it is the part that is ultimately responsible for the radiative damping. 

Unlike the dissipative long range part, the dominant contribution considered here 
is strictly conservative: one of the main results of this work is the demonstration that 
for a generic (not just a transonic) elastic string model the divergent self gravitational 
stress energy contribution given ("on shell") by the new formula]^] is derivable from 
a corresponding ("off shell") action contribution that is precisely the same as what is 
obtained as the four dimensional specialisation of a more general (higher dimensional) 
action formula that has recently been derived using an entirely different approach by 
Buonanno and Damour [jl^. The net effect is thus describable as an action renormalisation, 
whose effect is trivial in the case a Nambu-Goto model for which (consistently with what 
was suggested by exact analytic considerations in static configurations |]l8[) it simply 
vanishes, but non-trivial in the more general string models considered here. 



2 Gravitational force density 

Since the cosmic strings to which this work applies will be characterised by a gravitational 
coupling constant that is very small, Gm^ < 10^^ (where G is Newtons constant) even for 
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the heavyweight case of GUT strings, it is sufficient to base the analysis on hnearised grav- 
itational theory, as expressed in terms of a small but perhaps rapidly varying perturbation 
hpLv = ^Qfiu of a slowly varying cosmological background metric g^j,y . Moreover for the pur- 
pose of analysing the short range self interaction it will sufficient to neglect the background 
curvature altogether, i.e. to take unperturbed 4-dimensional spacetime metric Qfj^u to be 
flat, so that subject to the usual |jl^ DeDonder gauge condition, V^h^^y = ^Vyh^^, the 
linearised Einstein equations for the Eulerian metric perturbation will reduce to the well 
known form 

V^V'W = -8^TG{2f^"' - ffg''^) , (4) 

where T'^'^ is the stress momentum energy density tensor of the source. 

The problem of ultraviolet divergences for point particle or string models arises be- 
cause in these cases the relevant source densities are not regular functions but Dirac type 
distributions that vanish outside the relevant one or two dimensional world sheets. In 
the case of a string with local worldsheet imbedding given by = x^{a} in terms of 
intrinsic coordinates a* (i = 0, 1), so that the induced surface metric will have the form 
lij = dfiux'^,! x^,ji the relevant source distribution will be expressible using the terminology 
of Dirac delta "functions" in the form 

f^'^ = J T^'^ - x{a}] (fa , (5) 

where I7I is the determinant of the induced metric, and the surface stress energy density 
T^'^, is a regular tensorial function on the worldsheet (but undefined off it). 

In the simple elastic models considered here, the only internal field on the string will 
be a surface current density = e^-'x^^jVj) where ^ is a scalar stream function on the 
world sheet (which will be a free variable in the variation formulation described below), 
using the notation e*-' = —e^^ for the antisymmetric worldsheet measure tensor that is 
specified (modulo a choice of sign representing an orientation convention) as the square 
root of the induced metric, i.e. s^^Skj = 7\-. If there were a non-zero charge coupling 
constant q, as supposed in Witten's theory [|| of superconducting strings, then this would 
correspond to an electric surface current density = qc^. However the present discussion 
is concerned just with the early regime in which effects of electromagnetic coupling are 
negligible compared with those of gravitation. This means that in the relevant action for 
the "bare" string model (i.e. the non self interacting limit), as given by an integral of 
the form, I = J C ||7||^''^ (Pa, the specification of the relevant Lagrangian scalar C on the 
worldsheet will be given by a master function. A, that depends only on the undifferentiated 
background metric and the gradient of the stream function ■0, according the formula 

Z = A + iT^^/i^, , (6) 

in which the coefficient for the linearised gravitational adjustment term here is the surface 
stress-energy tensor that specifies the gravitational source in (^, which is given by T^'^ = 
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Since the ensuing field equations will evidently involve gradients of the stress energy 
tensor, their formulation will require the introduction of the appropriately defined hyper- 
Cauchy tensor (a relativistic generalisation of the Cauchy elasticity tensor of classical 
mechanics) which is defined by C^^^'^ = \\'^\\^^/'^d{TP'^\\'^\\^/'^) / dg^^p. In terms of this 
quantity, the dynamical equations obtained from the Lagrangian (^) can be shown ||l^ to 
be expressible in the standard form 

where the effective gravitational force density vector is given by 

using the notation = rj^Vy for the tangentially projected gradient operator, where r]^ 
is the tangential projection tensor, i.e. the index lowered form of the first fundamental 
tensor of the worldsheet, which is obtained simply by mapping its internal metric onto the 
spacetime background according to the formula rj^^'^ = ^^x'^^i x'^j . 



3 Renormalisation 

The problem with the application of (^) is of course that the linearised gravitational field 
will consist not just of a well behaved long range contribution, hfj^i^ say, but also of a 
divergent short range self interaction contribution, hfj_i, = h^y — h^^, that needs to be ap- 
propriately regularised in the manner recently described in the analogous electromagnetic 
case This routine procedure leads to a result that is proportional to the relevant 

source in (|^, which gives 

h^^ = 2GT{2%^-%''g^^), (9) 

where, as usual for a string self interaction in 4 dimensions, the proportionality factor has 
the form 

r=ln{AV5*n (10) 

in terms of an "ultraviolet" cut off lengthscale 5* representing the effective thickness of 
the string, and a much larger "infrared" cut off A, given by a lengthscale characterising 
the large scale geometry of the string configuration. As pointed out in the electromagnetic 
case 1 12 1, the corresponding regularised value of the gradient of such a divergent self field 
will be obtainable from the regularised self field by the application of the regularised 
gradient operator defined by 

% = %+^^^^., (11) 

where is the worldsheet curvature vector that is obtainable as the surface divergence 
of the fundamental (tangential projection) tensor, K^^ = Vj^ry^. 
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When the corresponding divergent self force contribution is evaluated substituting @ 
and (^) in (^), the result turns out, rather remarkably, to be describable as a renormal- 
isation of the stress energy tensor, since one obtains |jl^ a regularised self force that is 
expressible as world sheet divergence 

/g'^ = -v,f/^ (12) 

in which the relevant stress momentum energy density contribution from the gravitational 
self interaction has the form 

fg'^^ = K'^T'''' - jhp^TP^ri'"' + hf.X'"''"' ■ (13) 

It is also to be observed that this "on shell" self gravitational stress energy contribution is 
obtainable from a corresponding "off shell" self gravitational action contribution given by 

Ag = \m^u = ^GlilT^jT^^" - T/T/) . (14) 

This provides a regularised treatment in which the original worldsheet Lagrangian C is 
replaced by a regularised Lagrangian C = K-\- i^T^^h^i, involving only the well behaved 
part hfj^y of the gravitational field, whereby the divergent part h^jj is absorbed into a 
renormalised master function given by 

A = A + Ag . (15) 

The consistency of this treatment has been neatly confirmed by an independent investi- 
gation in which, by working entirely at the level of the "off shell" action in a space-time 



of arbitrary dimension, Buonanno and Damour |17] have recently obtained a general self- 
interaction formula whose specialisation to the case of gravitation in 4-dimensions is in 
precise agreement with the result (^) obtained here. 

For any simple elastic string model of the kind considered here, the master function A 
will depend just on the scalar magnitude that is specifiable as x = c^c^ = — 7*-^V',iV',j = 
—PfiP^ where the relevant momentum vector is defined by p'^ = V^ip = ^i'y'^^ip^j. Using 
the notation A' = dA/dx, it can be seen that the regularised self field will be given by 

h'"' = 4G r(2A'pV + X^'g'"" -AAl""), (16) 

using the notation -\J]j= g^y — r]^y for the (rank 2) worldsheet orthogonal projection tensor. 
The corresponding self gravitational action contribution (|l^ will be expressible as 

K^ = \G^{U -Tf = 2G^{xA'f . (17) 



6 



4 Effect on stress tensor and propagation speeds 

In analogy with the standard formula 

T^^ = Kt]^^ + 2A'p'^p^ , (18) 

for the "bare" surface stress energy tensor, the corresponding self gravitational stress 
energy tensor ( |l3|) will be obtainable in the form 

f/"" = Ag'n^"' + 2A^ p^'p'' . (19) 

If the current is timelike - as can be assumed without loss of generality in the wiggly 
Nambu-Goto string approximation - so that we have x ^ Oi then the energy density 
and tension will be given for the bare model hy U = —A and T = 2xA' — A, while the 
corresponding extrinsic (wiggle type) and longitudinal (sound type) perturbations speeds 
Cg and will be given Q by 

c2 = 1 - 2xA7A , c2 = 1 + 2xA7A' . (20) 

An elastic string state is describable as supersonic, transonic, or subsonic, according to 
whether the difference, 

c^-cl = -2x{ln{AA'})' , (21) 

is positive, zero, or negative. In terms of these quantities the corresponding renormalised 
energy density (for a given value of the current magnitude as specified by x) will be 
obtainable directly from the renormalised action ( p!5[ ) as 

U = -A = U-^GTUT{c-^-c^f. (22) 

However the evaluation of the corresponding renormalised tension is not so quite so simple: 
it works out to be 

f = T+jGTUT (c-i - c^f{l + 2c2) . (23) 



5 The case of the transonic wiggly string model 



In the exact Nambu-Goto case [16| the master function is just a constant, A = — m , where 
(in units with the Dirac Planck constant h set to unity) the parameter m is the relevant 
Kibble mass scale, which is a constant that can be expected to be of the same order of 
magnitude as the Higgs mass scale associated with the underlying symmetry breaking for 
a string of the ordinary "cosmic" kind, representing a vortex type defect of the vacuum. 
In this case one simply obtains A' = 0, so the divergent short range self gravitational 
contribution vanishes. 
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However for the purpose of a course grained description on larger scale, it is appropriate 
to use the transonic model ^ to represent the smoothed average over shortwavelength 
wiggles in the underlying microscopic Nambu-Goto model. For this transonic string model, 
the relevant master function has the form 



K = —m\J m? — X ■, (24) 

with non zero derivative 

A' = -mV2A. (25) 

This gives a non-zero self gravitional contribution attributable to interaction at interme- 
diate range, i.e. distances large compared with the wiggle wavelength in the underlying 
macroscopic model but small compared with the smoothing length on which the macro- 
scopic description is based. A noteworthy example of the application of this model is to 
the case of wiggles that are of purely thermal origin, for which the state parameter x will 
be given[^, 19| as a function of the relevant temperature B by the formula 



from which it can be seen that corresponding wiggle propagation speed will be given by 

It is evident that the special transonicity property 

= (28) 



of the "bare" model will not survive in the renormalised model, as obtained using (24) from 
([l5| ) and (|T7|), for which (still working just to linear order in the gravitational coupling) 
the difference between the squared propagation speeds is found to be given by 



GmHic-' - cj((l + clf + i(l - clf) . (29) 



The manifest positivity (since < < 1) of this result shows that the "dressed" gravita- 
tionally self interacting wiggly string model will be of supersonic type. 



6 Orders of magnitude. 

The regime of applicability of the foregoing analysis is of course subject to limitations. To 
start with, for the validity of the linearised gravitation equation (^) on which the entire 
analysis depends, the weak coupling condition 

GmH (30) 
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must be satisfied. This requirement is not very restrictive, because even the heaviest kind 
cosmic strings that are commonly considered in cosmological apphcations, namely those 
arising from G.U.T. symmetry breaking, are characterised by Gm^ « 10~^. For other 
kinds, such as those arising from electroweak symmetry breaking, the value of Gm^ will 
be smaller still. The smallness of Gm? will be partially counterbalanced by the fact that 
the regularisation factor / will be large compared with unity, but since, according to (|lO|), 
it arises as a logarithm, it can never be extremely large. In nearly all cases that are likely 
to arise in practice it will satisfy / ^ 10^, so the requirement ( pO[ ) will be satisfied by a 
large margin. 

Since the rigourous justification |lT| of the description of the macroscopic effect of the 
wiggles by the simple elastic model (^) depends on the supposition that self intersections 
of the string can be neglected (since if loop formation were important a more elaborate 
non-elastic model would be needed) the validity of the model as a precise representation 
is limited to the regime in which the effective energy density of the wiggles, as formally 
defined by 

U ~ T 

e = — — c^U -m'^ , (31) 
is small compared with the intrinsic energy density of the string, i.e. 

e < , (32) 
which is interpretable in the thermal case as meaning 

e < m . (33) 

Since this is equivalent to the restriction 

1 - « 1 (34) 

it can be seen that the difference ( p9| ) will be given approximately by the simple formula 

c2-c2~8Gre. (35) 

In an application of this kind, the magnitude, A say, typifying the wavelength of the 
wiggles over which the averaging is taken will provide an appropriate choice for the infra- 
red cut off, i.e. it will be natural to take A ~ A. Similarly the magnitude, a say, 
characterising the amplitude of the wiggles will provide the corresponding value for the 
ultra-violet cut off, which will thus be given by 5^, ~ a, provided that, as will usually 
be the case, the amplitude of the wiggles is not even smaller than the microscopic string 
radius, r say, in which case the latter would itself provide the relevant ultraviolet cut off. 
When the string is of the usual "cosmic" variety, representing an underlying vortex type 
defect of the vacuum, one expects the radius to be of the same order of magnitude as the 
Compton wavelength associated with the relevant Kibble mass, i.e. one expects to have 
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r ~ m^^. In most applications of interest, the relevant amplitude a will be considerably 
larger than this, and even for wiggles of purely thermal origin it will never be smaller: 
for the wiggles produced by a given temperature 0, as given by (p6|), one can estimate 
the relevant magnitudes as A ~ 0^^ and a ~ m~^, i.e. independently of the temperature 
the relevant amplitude will be of the same order as the cosmic string radius, a ~ r. This 
means that, in the usual physical applications, it will generally be possible to take the 
relevant cut-off ratio to be given by A/(5* ^ X/a. 

Since the effective energy density, e, contributed by small wiggles of wavelength A and 
amplitude a can be roughly estimated as e ~ rr?o? j}? ^ we see that the squared cut-off 
ratio appearing in the logarithm in ( p!o| ) will be given by the corresponding rough estimate 

/^i ~ m^/e, and hence that a reasonably accurate description should be obtainable by 
taking the regularisation factor itself to be given by an approximation of the form 



r^ln{^}, (36) 



where (e) is a constant chosen as some suitably weighted mean value of the energy density 
e = m^x/2A in the string segment under consideration. The requirement ( |3^ ) that the 
wiggle energy density should be small compared with the intrinsic energy of the string 
implies that the ensuing factor I will be reasonably large compared with unity, and the 
consideration that the dependence is logarithmic means that the result will be inensitive 
to the details of the particular prescription chosen to specify (e) . In order to obtain higher 
accuracy , one might be tempted to replace the fixed mean value (e) in (|3^) by the variable 
local value of e, which, by (^^, would be equivalent to taking the self gravitational action 
adjustment to be 

Ag~2Ge2ln{!^} (37) 

with 

However the appearance of improvement provided by such use of a variable rather than 
constant value for the renormalisation factor I is rather illusory, since the preceeding 
demonstration of renormalisability - as embodied in the formulae (39), (40) and (67a) 
- was dependent on the postulate that / should be constant. Moreover it can easily be 
checked explicitly that the only effect of the use of the apparently more precise formula 
(|37|) on the final formula (^5|) will be to replace the factor ^ by / — |. This adjustment 
would be significant only if e were so large as to be comparable with m^, so that the 
accuracy of the treatment would in any case be affected by other complicating processess 
such as loop formation due to self intersections. 
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